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Abstract 

The electronic structure of VO2 is studied in the frameworks of local density 
approximation (LDA) and LDA+J7 to give a quantitative description of the 
metal-insulator (MI) transition in this system. It is found that, both struc- 
tural distortion and the local Coulomb interaction, play important roles in the 
transition. An optical gap, comparable to the experimental value has been 
obtained in the monoclinic structure by using the LDA+C7 method. Based on 
our results, we believe that both, the Peierls and the Mott-Hubbard mecha- 
nism, are essential for a description of the MI transition in this system. 
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The observation of a metal-insulator (MI) transition in VO2 has attracted much interest 
for several years. Actually, this material has been considered a prototype of transition-metal 
oxides with an interesting MI transition. The transition occurs at temperature about 340 
K with a cell-doubling phase transition from the rutile to the monoclinic structure, its con- 
ductivity jumping by more than five orders of magnitude [I]. The transition is associated 
with a small change in volume, but significant internal atomic rearrangement characterized 
as both a pairing and an off-axis displacement of alternate vanadium atoms. In its insu- 
lating phase, VO2 has no magnetic order. The discovery of this transition initiated several 
experimental [@-|5|] and theoretical investigations [|BHTl"ll, trying to explain the explanation of 
the mechanism of this phenomenon. 

It is commonly accepted that Peierls and Mott-Hubbard pictures are two main mecha- 
nisms for MI transitions in half-filled band systems. Since VO2 is of a half-filled band type, 
and possesses both a structural transition as well as a significant Coulomb interaction (U) 
in the 3d bands of vanadiums, the question arises which mechanism drives the transition. 
This problem has been debated for nearly four decades, and continues to be a controversy. 
An early attempt was given by Goodenough ||, who, due to the existence of V-V pairings 
in the monoclinic structure, argued the transition to be of Peierls type: in his point of view, 
the V-V pairing and the pronounced octahedral distortion are responsible for the insulating 
state, and the gap is solely determined by the crystalline distortion. Somewhat later, Zyl- 
bersztejn and Mott proposed another model 0. With the information obtained from the 
properties of Vi-xCt x 2 , they argued that the band gap should essentially be due to the 
Hubbard mechanism. Nevertheless, both of the two models are based on qualitative analysis 
and have not offered us a satisfied scenario of the MI transition. 

Many theoretical efforts have been made to give a quantitative description of the MI 
transition in VO2 from first-principles calculations. Earlier works include a semiempirical 
band study by Caruthers et al. || and a cluster calculation on V-0 octahedrons by Sommers 
et al. ||. Most noticeable, in 1994, an ab initio molecular dynamics (MD) study on VO2 



performed by Wentzcovitch et al. |10j gave a good description to the structural performance 



of VO2, in which V-V pair distances along the c-axis in the monoclinic phase are further 
reduced from the experimental value. However, their calculated optical gap of -0.04 eV 
contradicts the observed value of 0.7 eV 0. A recent study based on the augmented spherical 
wave (ASW) method within the atomic sphere approximation (ASA) gave nearly the same 
band picture of electronic structure [TT| . 



Thus up-to-date results from LDA and other calculations can not provide a satisfactory 
quantitative description to the MI transition in VO2. We suspect that this failure is mainly 
due to the neglect of the local Coulomb interaction, which may be significant in 3d transition 
oxides. Since no full potential, all-electron calculations, and no quantitative study consid- 
ering a finite Coulomb interaction U have been reported so far, the driving mechanism is 
still unclear, and a comprehensive study with taking into account both effects is in order. 
In this article, we propose to study the mechanism of the MI transition in VO2 by using the 



LDA as well as the LDA+?7 method |12| , [13| . The latter has been proved to be an efficient 
and reliable tool to study systems with strong Coulomb interactions. Our results show that 
both, lattice distortion, and Coulomb interaction play important roles in the transition. It 
seems that correct and quantitative description can not be obtained without taking into 
account either of them. 

Our calculations are based on density functional theory within LDA and LDA+£7. 
The LDA calculations are performed in the full-potential linearized augmented plane-wave 
method (FP-LAPW) with local orbital extensions ||14| , and the full-potential linear muffin- 
tin orbital approach (FP-LMTO) jnj, while LDA+[7 calculations are carried out within 
the framework of the FP-LMTO approach ||15|| . The non-overlapping muffin-tin spheres are 
adopted in both approaches. For the FP-LMTO approach, three re's basis sets, -0.1, -1.0, 
and -2.5 Ry, are chosen. The potential is expanded in spherical harmonics inside the sphere, 
and Fourier transformed in the interstitial region. Therefore, no shape approximations are 
made to the density or potential. In both methods, the calculated results are checked for 
convergence with respect to the number of k points and the plane wave cutoff energy. In 
the LDA+[7 method, the screened U parameter used are ranging from 1.2 eV to 5.0 eV, 



which have been suggested by different authors, and are reasonable in early transition oxides 
PJ5| JT^Jr7[1 . The screened J is estimated to be 0.12 U. 

Our calculations have been carried out in the rutile structure of VO2 (space group 
P^/mnm) as well as in the monoclinic Mi structure (space group P2i/c). The lattice 



parameters for these two structures used in our work are taken from McWhan et al. [I8| and 
Longo et al. [f[9|| , respectively. In contrast to previous band calculations, our calculations 
were performed with full-potential methods. 

The densities of states (DOS) of the two structures, as determined by the LDA calcu- 
lation, are shown in Fig. |l|. It can be seen from Fig. ||(a) that the DOS in the metallic 
phase near the Fermi level is contributed mainly from vanadium 3d states, which are split 
by the octahedral field into two groups, a lower triply degenerate state of t 2g , and an upper 
doubly degenerate state of e g . The degeneracy of t 2g (e g ) is further eliminated due to the 
orthorhombic field. Two distinct changes are obvious when comparing the two phases. One 
is an upward shift of so-called tt* bands (d xz and d yz ) above the Fermi level, the other is a 
splitting of the d\\ band {d xy ) of about 1.7 eV. We argue that the former change is mainly 
caused by the V-V pairing. According to our calculation, the shifts of d xz and d yz band 
are about 0.18 and 0.27 eV with respect to the rutile phase, respectively. The occupation 
number of the ir* band is reduced from 0.45 e/atom in the rutile phase, to 0.15 e/atom 
in the monoclinic phase. These two effects greatly eliminate band crossings, and lead to a 
reduction of the DOS at the Fermi level (see Fig. |l](b)). Our calculations are in qualitative 
agreement with Goodenough's model || and quantitatively confirm the recent results by 
Eyert [|TTJ. The results obtained from the ab initio MD gave a smaller overlap of bands near 
the Fermi level than ours: in the MD calculation, the potential adopted is a pseudo one and 
the parameter used is theoretically optimized, i.e. the c-axis V-V pairing has been reduced 
from the experimental value of 2.62 A to 2.58 A. 

Our calculations show, confirming a suggestion by Eyert [[□]], that further reducing 
the V-V pairing to ~ 2.50 A even produces a non-overlap band structure; which even 
though the value is unrealistic, it is an interesting observation in itself. Our FP-LAPW 



calculations also show that artificial reduction of V-V paring have raised the total energy of 
the monoclinic phase comparing to that of the experimental lattice parameter. Specifically, 
the total energy in the experimental parameter is calculated to be of ~ 6 mRy/molecule lower 



compared to the case when the optimized value is used P0| . We also used the generalized 
gradient approximation (GGA) to calculate the band structure of VO2, with essentially no 
improvement. We conclude that the structural distortion plays an important role in MI 
transition in VO2, and it strongly alters the band structure of the insulating phase near 
the Fermi level. The upward shift of the n* and the separation of the upper d\\ band from 
the lower d\\ band caused by the V-V pairing leaves the lower d\\ band nearly filled. Even 
so, it can only lead to a reduction of carriers (or a small gap) near the Fermi level. The 
electronic properties of V0 2 in the insulating state, therefore, can not be described by a 
LDA band calculation. A significant improvement is obtained, as we demonstrate below, by 
supplementing the LDA by a local Coulomb term U. 

In order to get a better picture of the MI transition in VO2, we go beyond the LDA 
and check the part of correlation played in this system. The band structure of VO2 in the 
insulating phase has been computed by the LDA+£7 method, and is shown in Fig. ^. [U =0 
represents the result for LDA.] We see that the band structure of the monoclinic phase is 
indeed modified by U. The Coulomb interaction U has two effects on the band structure, 
see, for example, Fig. 0(b): the first effect is to cause a dramatic shift of the ir* band from 
the lower du band for ~0.7 eV, the other one is to cause a further separation of the du 
bands. The former effect is responsible for opening a gap; the latter one is not so important 
since the separation of the d\\ band by the structural distortion has been large enough, i.e., 
greater than the optical gap of 0.7 eV. Notice that this can only be achieved on the basis 
of the structural distortion. To illustrate this point of view, we have performed a LDA+[7 
calculation on the metallic structure of VO2, with the result that a gap is not opened by 
a value of U up to 7 eV. In this case, the overlap of the tt* and the d\\ band has not 
been eliminated, and the du band near the Fermi level is still of a half-filled type. Further 
increasing U is not reasonable in this system. This clearly indicates that the metallic phase 



is determined by the structure, and can be adequately understood by LDA theory 

By varying U in the reasonable range from eV to 5 eV in the insulating phase of VO2, 
we obtain gap values from -0.17 eV to 0.96 eV. The band gap as a function of the Coulomb 
interaction U is displayed in Fig. [| [A negative gap represents an overlap of the d\\ and the 
71* bands near the Fermi level.] A gap is opened, and exhibits roughly a linear relationship 
with U in the range of U between 1.6 eV and 5 eV. Specially, for U=4 eV, our calculation 
gives a gap of 0.67 eV, which is in good agreement with the experimental value of 0.7 eV. 
It is instructive and necessary to apply the LDA+£7 method to M0O2 and check whether 
this method are valid in our study. Although M0O2 crystallizes in the same monoclinic 



structure as low temperature VO2, it is still metallic. The LDA+U band calculations |20 
shows that the lower du band is completely occupied, and that the Fermi level is located in 
the lower-middle part of the ir* band. This condition is quite different from low temperature 
VO2. Since the overlap of bands near the Fermi level is about 1 eV, even a rather large U 
(up to 6.0 eV) can not open a gap in this system. This explains well why M0O2 is a still 
metal in spite of large U values. 

Figure f| shows the calculated DOS by LDA+U compared with the ultraviolet photoe- 
mission spectroscopy (UPS) spectrum || in the insulating phase of V0 2 . The UPS spectrum 
clearly shows that the 2p and 3d band are in the energy range of -2 to -9 eV, and to -2 
eV, respectively. We find that there is an overall good agreement between experiment and 
theory, with the exception of an underestimation of the 3d downshift by about 0.9 eV. Also, 
the calculated DOS fits well with the x-ray photoemission spectroscopy (XPS) and UPS 
spectra by Goering et al. || and Bermudez et al. j| First, the calculated DOS drops to zero 
at the Fermi level, which well explains the structure of UPS spectrum near the Fermi level. 
Second, the three peaks in the DOS curve centered at about -3.5 eV, -5.2 eV, and -7 eV 
are in good agreement with the O 2p structure in the UPS spectra. We also compared our 
result with the O Is x-ray absorption spectroscopy (XAS) spectra |21| which reflects the 
unoccupied states. The three peaks of it*, du and cr* can also be well explained by our DOS 
results. 
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Further discussion about 3d transition oxides such us Ti02 and &O2 will help us to 
understand the MI transition in VCV It is interesting to note that the above oxides in 
rutile structure exhibit different properties, i.e., Ti02 is an insolator, VO2 undergoes a MI 
transition and &O2 is a semi-metal. A simple model derived from the band picture can give 
a good explanation. There are two MO2 (M=Ti, V, Cr) units per cell in each systems while 
the numbers of occupied valence electrons in the du bands are quite different. In Ti02, since 
the O 2p orbitals are filled, the zero number of d electrons leave unoccupied d bands above 
Fermi level and cause an insulating phase. This is quite different from the case of VO2. There 
are two half-filled bands near the Fermi level which can not be split by a rather large value 
of U. This situation can be changed by cell-doubled distortion during the MI transition. In 
the insulating phase there are four VO2 molecular per cell which nearly fills the lower du 
bands. According to this picture, one may expect &O2 to be an insulator since the two 
bands near the Fermi level are fully occupied by four electrons, and the Coulomb interaction 
will separate these bands from other d bands. However, the ferromagnetic mechanism will 
account for the semi-metal property P^P^fl . 

In conclusion, the band structure of VO2 is calculated with the LDA and LDA+£7 meth- 
ods. Our results show that both, lattice distortion and Coulomb interaction, play important 
roles in the MI transition in VO2. First, the distortion strongly alters the band structure 
of the low temperature monoclinic phase relative to that of the high temperature rutile 
phase. It induces a small overlap of bands at the Fermi level but still leaves the system to 
be metallic. Second, the Coulomb interaction further splits the n* bands from Fermi level 
via the Mott-Hubbard mechanism. This cooperative effect causes an optical gap and thus 
the MI transition. With a reasonable value of U ~ 4 eV, we obtain a gap of 0.67 eV, which 
is close to that observed in experiment. A similar calculation is also performed in the rutile 
phase. There, no gap is found for U up to 7 eV. In conclusion, neither the Peierls nor the 
Mott-Hubbard mechanism provides a comprehensive picture of the transition. Our calcula- 
tion indicates that it is rather a combination of both, which gives a quantitative description 
of the metal- insulator transition in VO2. 
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FIGURES 

FIG. 1. Partial V 3d densities of states (DOS) of VO2 in (a) the high temperature rutile and 
(b) the low temperature monoclinic structure, respectively, within the LDA method. 

FIG. 2. Band structure of monoclinic V0 2 from (a) LDA and (b) LDA+C7, with U=4.0 eV. 

FIG. 3. Band gap vs screened (local) Coulomb interaction. A negative "gap" means that an 
overlap of the d\\ and the ir* bands. 

FIG. 4. Comparision of ultraviolet photoemission data [2] with the calculated DOS of mono- 
clinic (insulating) VO2. 
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(a) rutile V02 
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(b) monoclinic V02 
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